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Abstract 

A subcritical branching process in random environment (BPRE) is considered whose 
associated random walk does not satisfy the Cramer condition. The asymptotics for 
the survival probability of the process is investigated, and a Yaglom type conditional 
limit theorem is proved for the number of particles up to moment n given survival to 
this moment. Contrary to other types of subcritical BPRE, the limiting distribution is 
not discrete. We also show that the process survives for a long time owing to a single 
big jump of the associate random walk accompanied by a population explosion at the 
beginning of the process. 
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1. Introduction and Main Results 

In this paper we consider the asymptotic behavior of a type of subcritical branching 
processes in random environment. More specifically, the random environment is given 
by a sequence of independent and identically distributed (i.i.d.) probability distribu- 
tions on nonnegative integers, denoted by 7r = {ii n , n > 0} where 

oo 

7r„ = KViVi 2 \...},7r«>0, £>« = !, 
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which are defined on a common probability space (ft, A, P). Moreover, for a given 
environment it = {ir n }, the branching process {Z n , n > 0} satisfies 

Z„ = l, E w [s z " +1 | Zo,^,...,^] = (/ n ( s )) z ", 

where f n (s) = Yl < iLo' K n > 31 * s the generating function of 7r n ; in other words, ir n is the 
(common) offspring distribution for the particles at generation n. Here and below we 
use the subscript it to indicate that the expectation (or probability with the nota- 
tion P^) is taken under the given environment tt. As is shown in various articles on 
the branching processes in random environment, the asymptotic behavior of {Z n } is 
crucially affected by the so-called associated random walk {S n } defined as follows: 

S = 0, S n := X 1 + X 2 + • • • + X n , n > 1, 

where 

X n = log/^ 1 (l),n=l,2,.... 
are the logarithmic mean offspring numbers. For notational ease let 

f(s) = f (s), and X = X 1 . 

We call {Z n } a subcritical branching process in random environment if 

E[X]=:-a<0. (1) 

By the SLLN this implies that the associated random walk {S n } diverges to —00 almost 
surely, which, in view of the inequality P„. [Z n > 0) < [Z n ] = e Sn , leads to almost 
sure extinction of {Z n } . 

Subcritical branching processes in i.i.d. random environment have been consid- 
ered in a number of articles, see, for instance, Afanasyev (1980), Dekking (1988), 
Liu (1996), Souza and Hambly (1997), Afanasyev (1998), Fleischmann and Vatutin 
(1999), Afanasyev (2001), Guivarc'h and Liu (2001), Geiger et al. (2003), Vatutin 
(2003), Afanasyev et al. (2005), Bansaye (2008), Bansaye (2009), Afanasyev et al. 
(2010), and Afanasyev et al. (2011). According to these papers, a subcritical branch- 
ing process in random environment is called weakly subcritical if there exists 6 G (0, 1) 
such that E [Xe ex ] = 0; intermediately subcritical if E [Ae x ] = 0; and strongly 
subcritical if E [Xe x ] < 0. 

The classification above is not exhaustive, though. An important exception is that 
the random variable X is such that E [Xe 6 ^] = 00 for any 9 > 0, and this is the focus 
of the present paper. To be more specific, we suppose that a 2 := Var(X) < 00 and, in 
addition, following the custom of writing / ~ g to mean that the ratio j jg converges 
to 1, we have, as x — > 00, 

A(x) := P (X > x) ~ for some p > 2, (2) 
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where l(x) is a function slowly varying at infinity. Thus, the random variable X does 
not satisfy the Cramer condition. 

We will see below that for the non-Cramer case, similarly to other cases, the asymp- 
totics of the survival probability and the growth of the population size given survival 
are specified in main by the behavior of the associated random walk. However, the in- 
fluence of the associated random walk for the non-Cramer case has essentially different 
nature: for the Cramer cases, the survival for a long time happens due to the "atypi- 
cal" behavior of the whole trajectory of the associated random walk that results in its 
smaller, then usually, slope for the strongly subcritical case (Afanasyev et al. (2005)), 
in its convergence to a Levy process attaining its minimal value at the end of the 
observation interval for the intermediately subcritical case (Afanasyev et al. (2011)), 
and in the positivity of its essential part for the weakly subcritical case (Afanasyev 
et al. (2010)). For the non-Cramer case, the process survives for a long time owing 
to a single big jump of the associated random walk at the beginning of the evolution 
which, in turn, is accompanied by an explosion of the population size at this moment; 
see Lemmas 14 and 15 for the precise information. Besides, the number of particles at 
a distant moment n given its survival up to this moment tends to infinity for the non- 
Cramer case, while for the other types of subcritical processes in random environment 
such conditioning leads to discrete limiting distributions with no atoms at infinity. 

One of our assumptions is the following (technical) condition for A(x): for any fixed 



Next, for any offspring distribution n = {n^ : i > 0} with generating function f(s), 
denote 



h > 0, 




(3) 



x 





Introduce the following 



Assumption 1. (i) There exists 5 > such that, as x — > oo, 




(ii) As x —> oo, (under probability P,) 



e~ x ) \X > x 




(5) 



where 7 is a random variable which is less than 1 with a positive probability. 
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It can be shown that if ttq is either almost surely a Poisson distribution or almost 
surely a geometric distribution, and if (2) is satisfied, then (5) holds with 7 = 0. More- 
over, it is not difficult to give an example of branching processes in random environment 
where 7 is either positive and less than 1 with probability 1, or random with support not 
concentrated at 1. Indeed, let 7 be a random variable taking values in [0, 1 — 5] C [0, 1] 
for some 5 G (0, 1], and p and q,p + q — l,pq > 0, be random variables independent 
of 7 such that the random variable X := log(l — 7) + \og(p/q) meets conditions (1) 
and (2). Define f(s) := 7 + (1 -7)9/(1 - ps). Then /' (1) = (1 ~i)p/q = exp(X), 
and it is straightforward to show that for any e G (0, 1), 

lim P (|/(1 - e~ x ) -i\>e\X>x) < lim P (X - x < - loge | X > x) = 0, 
therefore (5) holds. 

Let us briefly explain Assumption l(ii). For any fixed environment 7r and any x > 0, 
let C n (Zie~ x ) be the distribution of Z\e~ x . Note that this actually only depends on 
7r . We will show in Lemma 9 that if condition (2) holds, then (5) is equivalent to 
the following assumption, concerning weak convergence of random measures: (under 
probability P,) 

conditional on {X > x}, C n (Zie~' x ) =^> j5 + (1 — 7) 5^ as x — > 00, 

where 5q and 5^ are measures assigning unit masses to the corresponding points. 

In what follows we assume that the distribution of X is nonlattice. The case when 
the distribution of X is lattice needs natural changes related to the local limit theorem 
that we use in our proofs (see Proposition 5 below). 

Define 

fk,n(s) := /fc(/fc+i(...(/„_i (s))...))> < k < n- 1, and /„,„ (s) := s. (6) 
When k — 0, /o, n (s) = E 7T (s Zn ) is the conditional probability generating function of 

The following is our first main theorem which deals with the the survival probability 
of the process. 

Theorem 2. Assume conditions (1), (2), (3) and Assumption 1. Then the survival 
probability of the process {Z n } has, as n — >■ oo ; the asymptotic representation 

P (Z n > 0) ~ KP (X > no) , 

where 

00 

K:=^E[l-/ 0j ( 7 )]G(0,oo), (7) 

3=0 

and 7 is a random variable that has the same distribution as the 7 in Assumption l(ii) 
and is independent of the underlying environment {n n } (and consequently o/{/oj} ). 
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As we mentioned earlier, if ttq is either almost surely a Poisson distribution or almost 
surely a geometric distribution, then (2) implies (5) with 7 = 0, so the constant K 
becomes 5^1 P(Zj > 0). In this case we can give the following intuitive explanation 
of Theorem 2: Let 

U n = inf {j : X 3 > na} (8) 

be the first time when the increment of the random walk S := {Sj, j > 0} exceeds na. 
Then one can show that the event {Z n > 0} is asymptotically equivalent to {U n < 
n, Zu n -\ > 0} = Uj <n {Zj-i > 0, U n = j}. Now for each fixed j > 1, P(Zj_i > 0, U n — 
j) ~ P(Zj_i > 0)P(X > na), and hence, not rigorously, 

00 

P{Z n > 0) ~ P{U n < n, Z Un ^ > 0) ~ P ( Z i-i > °) ' P ( X > na ) = KP{X > na). 

3=1 

In fact, we may say more: (even in the general case when 7 ^ 0,) the process survives 
owing to one big jump of the associated random walk which happens at the very 
beginning of the evolution of the process; moreover, the big jump is accompanied by 
a population explosion which leads to survival. See Lemmas 14 and 15 for the precise 
information. 

The next result gives a Yaglom type conditional limit theorem for the number of 
particles up to moment n given survival of the process to this moment. Recall that 
a 2 = Var(X), and U n is defined in (8). 

Theorem 3. Assume conditions (1), (2), (3) and Assumption 1. Then for any j > 1, 
lim P(U n = j I Z n > 0) = E(l - foj-^/K. 

n— s-oo 

Moreover, 

£ (y — < t < 1 z n > oW (1, < t < 1) , (9) 

V Z Un exp{b [nt]vUn - b Un ) J 



C {^j= (log (Z [nt]vU jZ Un ) + nta) , < t < 1 
and for any e > 0, 

£ (~/^ ( lo S ( Z [nt]/Z [n£] ) + n(t - e)a) , e < t < 1 



Z n > ) (B t , < t < 1), 



Z n >0) =>> (St-S e , e<t<l), 



where the symbol means weak convergence in the space D[0, 1] or D[e, 1] endowed 
with Skorokhod topology, and B t is a standard Brownian motion. 

Therefore after the population explosion at time U n , the population drops exponen- 
tially at rate a, with a fluctuation of order exp(0(\/k)) with k the number of generations 
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elapsed after the explosion. Moreover, it follows from (9) and the continuous mapping 
theorem that 

C (log (Z [nt]vUn /Z Un ) - (S [nt]vUn - S Un ), < t < 1| Z n > 0) =^ (0, < t < 1), 

and, therefore, after the big jump, at the logarithmic level the fluctuations of the 
population are completely described by the fluctuations of the associated random walk. 

2. Some Preliminary Results 

We list some known results for the random walk S and establish some new ones. 
Define 

M n = max Sk L n = min Sk, r n = min (0 < k < n : Sk — L n \ , 

l<fc<n ' 0<fe<n 

t(x) = inf {k > : S k < -x} , x > 0, 
and t = r(0) = inf {k > : S k < 0}. Further, let 

k=i 

which is clearly finite given conditions (1) and (2). 

Proposition 4. [Borovkov and Borovkov (2008, Theorems 8.2.4, P a 9 e 376)} Under 
conditions (1) and (2), as n —> oo, 

P (L n > 0) = P (r > n) ~ e D P (X > an) . 

Next, let Y = X + a. Then Y is a random variable with nonlattice distribution, 
and with zero mean and finite variance. Moreover, as x — > oo, the function 

B{x) := P(Y > x) (= P(X > x - a) = A(x - a)) ~ /3 > 2, 

and satisfies a modified version of (3) by replacing A(x) with B(x). 

Proposition 5. [Borovkov and Borovkov (2008, Theorem 4-7-1, page 218)} Assume 

(2) and (3). Then with S n := Y 1 + ... + Y n , where Y i = Y and independent, we have 
for any h > 0, uniformly in x > N \/n\og(n + 1), as N — >■ oo ; 

p(S.e[«,» + fc)) = ^W(i + o(l)). 
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The uniformity of o(l) above is understood as that there exists a function 5(N) J, 
as N — > oo such that the term o(l) could be replaced by a function 8 h (x,n) with 
\S h (x,n)\<5(N). 

Based on Propositions 4 and 5 we prove the following 
Lemma 6. Assume conditions (1), (2) and (3). Then, as n — >■ oo ; 



E [e s "; r„ = nl = E [e 5 "; M„ < Ol ~ — P (X > an) , 

Tl 



where 



Ki := - exp J V -E [e 5 "; 5 n < Ol 1 < oo. (10) 
a \tl n J 

Proof. The first equality follows from duality. More specifically, the random walks 
: fc = 0, 1, . . . , n} and {S' k := — S n -k '■ k — 0, 1, . . . , n) have the same law, and 
the event {r„ = n} for {S*.} corresponds to the event {M' n < 0} for {S' k }. 

Next we evaluate the quantity 

E [e Sn ; S n < 0] = E [e 5 ™; - (£ + 2) logn < S n < 0] + O (n"^ 2 ) . (11) 
Clearly, for any h > 0, 

e -(fc+i)h . P (^_( k + ^ + an < 5 n < _£;/i + a n) 

0<k<(/3+2)h- 1 logn 

< E [e Sn ; — (13 + 2) log n < S n < 0] 

< £ e -^.p ^_(jfe + + an < S„ < -Jfe/i + an) . 



0<fc<(/3+2)/i- 1 logn 

By Proposition 5, in the range of k under consideration, as n — > oo, 

hpn 

(—(k + l)h + an) 
hj3_ 
a 

where o(l) is uniform in < k < ((5 + 2) n -1 logn. Now passing to the limit as n — >■ oo 
we get 



(-(k+ l)h + an < S n < -kh + an) = . - ^ -B (—(k + l)h + an) (1 + o(l)) 

V / ( — (fc + lj/i + an) 

A (an) (l + o(l)), 



< Ita irf »E[^ ; -(ff + 2)lcgn<S,<0] 

n^oo pA (an) 

^ aE[e s ";-(/3 + 2)logn<5 n <0] 

< hm sup — — — r - 

n^oo pA (an) 

oo 



< h J2 e 



kh 

e 
k=0 



Letting now h — > 0+ we see that 



Um aE [e s ";-(^ + 2)logn<g w <0] = i 
n-^oo (3 A (an) 

Combining this with (11) we conclude that, as n — > oo, 



E \e Sn ; S n < 0l = —A (an) (1 + o(l)) ~ -P(X > an). (12) 
L J a a 

Furthermore, we know by a Baxter identity that 



exp 



J2 -E [e 5 "; S n < oH = 1 + £>E [e 5 "; M n < 0] , 



n=l ) n=l 



see for example Chapter XVIII. 3 in Feller (1966) or Chapter 8.9 in Bingham et al. 
(1987). From (12) and Theorem 1 in Chover et al. (1973) we get 

E \e Sn -M n < Ol ~ — P(X > an), 

where K x is given by (10). That K\ < oo follows from (12). □ 
Corollary 7. Under the conditions of Lemma 6, the constant K in (7) is finite. 

Proof. Clearly, 

1 - /oj(7) < 1 " /oj(0) = P«(Zj > 0) = mm. P^Z, > 0) < min. e 5 » = e s ^ . (13) 
Thus 

oo oo 

K = ^E[l-/ ,( 7 )]<^E[e^] 

3=0 3=0 

oo j oo oo 

= E E E [ e5< ; T ^ = j ] = E E [ e5l 5 r ^ = j ] • E p (Vi > °)- 

j=0 i=0 i=0 j=i 

The last term is finite by Lemma 6 and Proposition 4. □ 

Next, recall that U n = inf {j : X, > na}, and r = infjj > : Sj < 0}. The next 
result says that if the associated random walk remains nonnegative for a long time, 
then there must be a big jump at the beginning. 

Proposition 8. [Durrett (1980, Theorem 3.2, page 283)] If conditions (1) and (2) 
hold then . 

lim P (U n = j\r >n) = — P (r > j - 1) . 

n— >oo rL/T 
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3. Proof of Theorem 2 

We first introduce the following convergence statements, which will be shown to be 
equivalent to each other: (under probability P,) as x — > oo, 

(i) for any function S(x) satisfying lim^oo S(x) — 0, C (/ (1 — exp(— x(l + S(x)))) \X > x) 

(ii) for any function S(x) satisfying lim^oo S(x) — 0, C (/ (exp(— exp(— x(l + S(x))))) \X > x) 

=► ^(7); 

(iii) for any A > 0, C (/ (exp (— A exp(— x))) \X > x) =>- £(7); and 

(iv) conditional on {X > x}, C n (Zie~ x ) =>- 7<5 + (1 — 7) 5^. 

Lemma 9. Assume condition (2). Then the convergences (i) ~(iv) above are equiva- 
lent, and are all equivalent to (5). 

Proof. We will show that (5)=^(i)=^(ii)=^(iii)«=> (iv), and finally (iii) => (5) . 

We first prove that (5) implies (i). By (2), the events {X > x} and {X > x(l + 
S(x))} are asymptotically equivalent to each other (in the sense that lim :c ^ 00 P(X > 
x\X > x(l + 5(x)) = \im x ^, OQ P(X > x(l + S(x))\X > x) = 1), hence it follows, for 
example, from Lemma 17 in Lalley and Zheng (2011), that 

lim P (f(l - e - x(1+5(x)) ) > y\X > x) = lim P (f(l - e -*(i+*(*))) > y \X > x(l + S(x))) 

= lim P (/ (1 - e~ x ) > y\X > x) . 

X-¥00 \ \ / 1 

Next we prove that (i) implies (ii). It is easy to see that for any 5(x) — > 0, for all 
sufficiently large x, 

1 - exp(-x(l + 5(x) + l/x)) > exp(- exp(-x(l + 5(x)))) > 1 - exp(-x(l + 8(x))), 
and consequently, by the monotonicity of /, 

f(l-exp(-x(l+8{x)+l/x))) > /(exp(-exp(-x(l+«J(a;))))) > f(l-exp(-x{l+6(x)))). 

Now, by (i), taking 5(x) to be 5(x) + l/x and 8{x) respectively, we have that both the 
first and the third random variables, conditional on {X > x}, converge in law to 7. It 
follows that the middle random variable also converges, implying (ii). 

To show (iii) from (ii) we simply take S(x) = —\og(X)/x. Moreover, (iv) and (iii) 
are equivalent since / (exp (— \e~ x )) is the Laplace transform of Z x e~ x . 

Finally we derive (5) from (iii). In fact, for all x sufficiently large, 
exp(— exp(— x)) > 1 — exp(— x) > exp(— e • exp(— x)). 
Hence, again by the monotonicity of /, 

/(exp(-exp(-a;))) > f(l - exp(-x)) > /(exp(-e • exp(-x))). 
The convergence in (5) then follows from (iii) by taking A to be 1 and e. □ 
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Corollary 10. Assume conditions (2) and (5). Then for any function 5 (x) satisfying 
lim-r-Hx, 5(x) = 0, 

E[ 7 ] = lim E [/(l - e- x{l+5{x ^))\X > x] = lim E [E^ \ e - XZl/e l \ X > x] . (14) 
In particular, 

E [7] = lim P(Zi < e x ( 1+ '( x ))|X > x). (15) 

a;— >oo 

Proof. This follows by applying the dominated convergence theorem to the conver- 
gences in (i) and (iii). □ 

Lemma 11. If (1), (2) and Assumption 1 are valid, then 

lim P (exp(— na — 2n/logn) < P w (Z n > 0) < exp(-na + n 2/3 )) = I. (16) 

In particular, for any sequence 5 n such that n(5 n — 2/ logra) — > 00, 

lim P (Z n > I Z > e n{a+Sn) ) = I. (17) 

n->oo 

Proof. The second claim (17) follows directly from the first one, so we shall only 
prove (16). We have (see, for instance, Geiger and Kersting (2000)) that 



n-1 '"• " 1 



V^Z n > 0) = ( e~ s " + Y,9k{fk + i,n^)y~ Sk J , (18) 

where 



fc=0 



< \ 1 1 

9k{s) :-- 



l-fk(s) f' k (l)(l-s) 
meets the estimates 

< g k (s) < 2rj k+1 with rj k+1 := rj(n k ). 
Introduce the events 

G n := j max |,S fe + A;a| < n 2 / 3 j and H n := j 1 + ^ r] k < 2ne n/log "| . (19) 

By the functional central limit theorem, lim n ^oo P (G n ) = 1. Further, by Assump- 
tion l(i), 

1 - P (H n ) < nP U > e n/logn ) = n x o I ^ ^ ) = o ( -^—) . 

V ; " V 7 I n(log(n/logn)) 1+<5 / \log'n/ 
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Thus, Hindoo P (H n ) = 1, and, consequently, 

lim P (G n DH n ) = l. (20) 

n— >oo 

It then suffices to show that on the event G n fl H n , 

exp(— na — 2n/ logn) < P n {Z n > 0) < exp(— na + n 2 ^ 3 ) 

for all sufficiently large n. 

The upper bound follows from the evident estimates 

PAZn > 0) < [Z n ] = cxp(S n ) < exp(-na + n 2/3 ). 

As to the lower bound, since gk(s) < 2n fc+1 , we have 

PjZ n > 0) > —i . 

e " + 2 Lfc=o Vk+ie bk 

Finally observe that on the event G n fl H n we have 

n— 1 / n— 1 

na+n 2 / 3 



e- 5 " + 2 J2 Vk + ie- % < 2e na+n2/3 1 + ]T 



fc+i 

fc=0 \ fc=0 

< 4 e "a+« 2/3 . ne ra / lo S™ < e na+2ri/ log n ^]_) 

for all sufficiently large n. □ 
Lemma 12. Assume (1), (2) and Assumption 1. Tnen /or any fceN, 
lim P (Z n > | Xi > na, Z = fc) = E [l - 7 fc l . 

n— >oo 

Proof. We only prove for the case when /c = 1. We have 

P (Z n > | X 1 > na) = E [P n (Z n > 0) | X 1 > na] = E [1 - / (/i,„(0)) | X x > na] . 
Write 

1 - / 1>n (0) = P v (Z n > Ofa = 1) = e—a+CW). 

According to the previous lemma, there exists a deterministic function 5(n) — > as 
n — > oo such that 

|C(n)|<5(n) 

with probability approaching 1 as n — >■ oo. The conclusion then follows from Corol- 
lary 10. □ 

Lemma 13. Assume conditions (1), (2) and (3). Then for any e > there exists M 
such that for all n sufficiently large, 

P (Z n > 0; r n > M) = E [P n (Z n > 0); r n > M] < eP (X > na) . 

11 



Proof. Using the inequality P n (Z n > 0) < e Srn established in (13), we have by Propo- 
sition 4 and Lemma 6 that 

n 

E [P„(Z n > 0); r n > M] < E [e 5 -; r„ = k] 

k=M 

= V[e Sk ;r k = k]P(L r ^ k >0)+ ^ E [e Sk ; r k = k] P (L n _ fc > 0) 

M<k<n/2 n/2<k<n 

< P (L [(n+1)/2] > 0) £ E [e s - r fc = fc] + C P {X > H ^ P (L fc > 0) 

fc=M " fc<n/2 

< eP (X > an) . 

□ 

Now we are ready to prove Theorem 2. 

Proof of Theorem 2. For fixed M and N we write 

P(Z n >0) = P (Z n > 0; U n < NM) 

+P (Z n >0;U n > NM; r n > M) + P (Z n > 0; U n > NM; r n < M) . 

By Lemma 13, for any e > there exists M such that for all sufficiently large n, 

P (Z n >0;U n > NM; r n > M) < P (Z n > 0; r n > M) < eP (X > na) . 

Moreover, by Propositions 8 and 4 there exists N such that for all sufficiently large n, 

P[Z„ :-():/; :- XM:t„ < M) < P (U n > NM; r n < M) 

M 

= Ys P ( U n>NM;r n = k) 

k=0 

M 

< ^2 P i T k = k)P (U n > (N — 1) M; t > n — M) 

k=0 

< (M + l)P(U n > (N- 1)M;t > n- M) 

< eP (r > n - M) < 2e D eP {X > na) . 

Hence the main contribution to P (Z n > 0) comes from P (Z n > 0; U n < NM), i.e., 
when there is a big jump of the associated random walk at the beginning. 

To proceed, we introduce the events 

A k = A k (n) := {Xi <na,l <i < k} , fc = l,2,...,n. (22) 
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For each fixed j, we have by the Markov property 

P(Z n > 0; U n = j) 
= P(Aj_i; Xj > na; Z n > 0) 

oo 

= P (A-i; = k ) P ( X i > na ) P {Zn-j+i > 0|Xi > na; Z — k) . (23) 
fc=i 

Clearly P(Aj_i, Zj_i = fc) = P(A,_i(n), = fc) increases to P(Zj_i = fc) as 
n — ?> oo. Dividing both sides of (23) by P (Xi > na) and applying the dominated 
convergence theorem and Lemma 12 yield 

P(Z >0;fl- n =j) = g P( ^._ i = A; . Vi) p ( ^ +i> o|x 1 >na,Z = A;) 

oo 

~ ^P(^_ 1 = A;)E[l- 7 fe ] 

k=l 

= E [1 — /oj-i (7)] ■ (24) 

The last equality holds because of the independence assumption that we put on 7 and 
{foj}- Hence 

00 

P (Z n > 0) ~ E t 1 " /o^-i (7)1 • P > H (25) 
as desired. □ 



4. Functional Limit Theorems Conditional on Non-extinction 

The proof of Theorem 2 shows that the process survives owing to one big jump 
of the associated random walk which happens at the very beginning of the evolution 
of the process. This motivates the study of the conditional distribution of U n given 
{Z n > 0}, which is the content of the next lemma. 

Lemma 14. For any j > 1, 

lim P(U n = j\Z n > 0) = E [1 - foj-M] /K, (26) 

n— >oo 

where K is as in Theorem 2. In particular, for any e > 0, there exists M > such 
that 

limsupP(£/ n > M\Z n > 0) < e. (27) 

n— >oo 

Proof. The first claim follows from the representation 

P{Un = jlZn>0) = ^M^Jl x ^>^ 



P(X x >na) P(Z n >0) 
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and relationships (24) and (25). 

Estimate (27) follows from (26) and Corollary 7. □ 

Thus, we have demonstrated that given survival to time n, there must be a big 
jump at the early time period. Next lemma complements this by showing that for 
survival of the process such a big jump will be accompanied by a population explosion. 

Let Zj(i) be the offspring size of the i-th particle existing in generation j — 1, and, 
as we shall deal with maxi<j<z [/n _ 1 Z Un {i) repeatedly, define 

Nu n -=, max Z Un (t). (28) 

Lemma 15. For any sequence h n such that h n < n and h n — > oo, and 5 n — > satisfying 
n(S n — 2/ logn) — >■ oo, 

lim P (U n < h n , N Un > e nia+s ^ | Z n > 0) = 1. 

Proof. We first estimate the probability 

P (U n < h n , N Un > e^^-Un > 0) 

oo 

= E P ^ = 3> Z i~l = k y P ( N Un > en(a+5n) I U n = J, Zj-! = k) 

j<h n k=l 

■ P(Z n > I N Un > e n( ~ a+5 -\ U n = j, Zj_, = k). (29) 
Recall the events Aj(n) that we defined in (22). As n — > oo, 

P{U n = j, Zj_! — k) — P(Z,-_i = k, Aj_i(n)) • P{X > an) ~ P(^_i = ife) • P(X > an). 
Moreover, by (a simple generalization of) (15), for any fixed j and k, 
lim P (N Un > e n ^ | U n = j, Zj.! = k) = E[l - 7 fc ]. 

Finally, by (17) in Lemma 11, 

P(Z n > | N Un > e < a+5 -\ U n = j, Z 3 . x = k) -+ 1. 

Dividing both sides of (29) by P[Z n > 0) and applying Theorem 2 and Fatou's lemma 
we get the conclusion. □ 

The arguments above lead to the following lemma, which says that conditioning on 
{Z n > 0} is asymptotically equivalent to conditioning on {U n < h n ,N Un > e n ( a +M}. 

Lemma 16. For any sequence h n such that h n < n and h n — > oo, and 5 n — > satisfying 
n(S n — 2/ logn) — > oo, 

||P(- \Z n >0)- P(- 1 U n < h n ,N Un > e< a+ ^)\\ TV ^ 0, (30) 

where \ \ ■ \ \tv denotes the total variation distance. 



14 



Proof. This follows from Lemma 17 in Lalley and Zheng (2011), relation (17), and the 
previous lemma. □ 

Lemma 16 shows that to prove the functional limit theorems for the population size 
up to moment n conditioned on survival of the process to this moment, we need only 
to establish limit theorems under the condition {U n < h n , Nu n > e n ( a +M}. 

Lemma 17. For all sufficiently large n, on the event G n fl H n (as defined in (19)) we 
have 

E "!fc\ < e Mna + 2n/ logn). 
Proof. Recall the generating functions /fc, n (-)'s that we defined in (6). Clearly, 

n-1 



fc=0 

and 

n-1 n-1 r I \\ n ~ l 

fo, n (s) = n mf»iM) • e r (f ! n m+iM)- 

i=0 k=0 Jk\Jk+i,nK & )) j=k+1 

Hence, letting s = 1 we get 

n-1 

E n [Z n (Z n - 1)] = /J B (1) = 2e 2S " -J2vk + ie- Sk+i . 

k=0 

Thus, 

n-1 

[Z 2 n ] = yj n (l) + [Z n ] = 2e 2S - -J2vk + ie- Sk+l + e s " 

k=0 

implying 

E. [Z 2 n ] 



721 n ~ 1 



exp(2S n ) 



2-J2vk+ie' Sk+1 +e-' s ". 



k=0 

The needed conclusion then follows from an argument similar to (21). □ 

For the following two lemmas we fix a sequence h n such that 

h n — > oo, and h n /n — > 0. (31) 

A simple example of such a choice is h n = logn. We also take 5 n to be a sequence 
satisfying 

S n ->■ and n(5 n - 2/ logn) ->■ oo. (32) 



15 



Lemma 18. Suppose (31) and (32) hold. Then 



Z[nt]VU n 



z u n exp(5 , [ nt ]v[/ n - S Un ) 



,0 < t < 1 



=>- (1,0 < t < 1). 



Proof. Let 7r' = {7T{/„ , Hu n +i, • • •} be the random environment after time U n , and S' m := 
Sm+u n — Su n be the random walk after time U n . Further, define G' n and H' n for the 
random environment ir' in the same way as in (19). Then, by (20), as n — > oo, G' n C\H' n 
occurs with probability approaching one, so we need only to prove the convergence on 
the event G' n H H' n . 

We first prove the marginal convergence, by a mean-variance calculation. Denote 
k = [nt\. By our assumption U n < h n for an h n /n — > 0, implying k = [nt] > U n for all 
t > and for all n big enough. Hence 

E n (Z k \Z Un ) = Z Un ■ exp(S' k _ Un ); 

moreover, by Lemma 17, for all n big enough, on the event G' n D H' n , 

Yai n (Z k \Z Un ) = Z Un ■ VaT w/ (Z k _ Un \Z = 1) < Z Un ■ exp(2S' k ^ Un ) ■ exp(na + 2n/ logn). 

Consequently, when N Un > exp(n(a + S n )) and, therefore, Z Un > exp(n(a + S n )), we 
have 



Var, 



^ exp(na + 2n/ logn) ^ exp(na + 2nj logn) 
Zi/ B exp(n(a + «J n )) 



Next, by Slutsky's theorem (see, e.g., Ferguson (1996)) we have convergence of finite 
dimensional distributions. Furthermore, since Z^ nt ^u n / {Zjj n exp(S , ' rat j v(7 _ v )) are mar- 
tingales (with respect to the post-C/ n sigma field 3~({ nt yj Un y where = a{-n; Zj,j <i)), 
to prove the convergence in the space D[0, 1] we need only to show, by Proposition 1.2 
in Aldous (1989), the uniform integrability of {Z [nt]vUn /(Z Un ew(S{ nt ]s/u n -uJ)} for an y 
fixed t. This follows from the above calculation, demonstrating that the elements of 
the martingale sequence in question are bounded in L? . □ 

The lemma just proved serves as an LLN type result; the following lemma gives the 
CLT type statement. 



Lemma 19. Suppose (31) and (32) hold. Then 

1 (log (Z [nt]vUn /Z Un ) + nta) , < t < 1 



U n <h n ,N Un >e n ^ a+s ^ 
=>(B u 0<t<l), 
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and for any e > 0, 

' U n <h n ,N Un >e n ^ 



C ( -i= (log (Z [nt] /Z [n£] ) + n(t - e)a) , e < t < 1 



=>(B t -B e ,e<t<l), 

where a is the standard deviation of X , and B t is a standard Brownian motion. 

Proof. By Lemma 18 and the continuous mapping theorem, 

£ (log (Z [nt]vUn /Z Un ) - (S [nt]vUn -S Un ),0<t<l\U n <h n , N Un > e n(a+5 ^) 

=>• (0,0 <t < 1). 

The conclusions then follow from the standard functional central limit theorem for the 
post-£4 random walk {S[ nt ] VUn - S Un }. □ 

Now we are ready to prove Theorem 3. 

Proof of Theorem 3. The first claim follows from Lemma 14. The second statement 
follows from Lemmas 16, 18 and 19. □ 
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